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1 Introduction and summary 

* (1) We study multi-flavored four dimensional QCD by adding soft supersymmetry breaking 

(SSB) terms to supersymmetric QCD (SQCD) 

* (2) Certain exotic phenomena discovered in SQCD persist also in the regime of small, 

compared to Aqcd, SSB parameters. 

* (3) The exotic behavior includes (i) a baryon number violating vacuum, (ii) massless 

composite fermions and (iii) Seiberg's duality 

* (4) In the decoupling limit most of the "exotic" phenomena disappear. 

* (5) We argue that all the exotic behavior of the region close to the SUSY limit should be 

detected in lattice simulations. 



2 On the effective action of SQCD 

SU{Nc) SQCD contains a gauge superfield Wa, a chiral quark and an anti quark superfields 
Qq, where i = 1, . . . ,Nf is a flavor index and a = 1, . . . , iV^i is a color index. These 
superfields transform under the quantum global symmetry 

SU{Nf)L X SU{Nf)n x U{1)b x U{1)r, (2.1) 

in the following way 
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The low energy effective action of SQCD is expressed, for small values of A^/, in terms of the 
meson superfield 

T,=Q'-Qy (2.2) 

Beginning at Nf = N^, Seff is built also from chiral baryon and anti-baryon superfields. The 
baryon 

R. — pai-'-a-Ncf . . . nh . . . nJNc (9 

transforms in the {Nf — Ai'c)-index antisymmetric tensor representation of SU{Nf)L, and, 
similarly, an antibaryon chiral superfield Hn^-Nc) i^^ji^ from Nc powers of the field Q. 

^We take here Nc > 3. The special case of A^c = 2 is discussed in ref.[^. 



Using the gauge supermultiplet, it is possible to build the (massive) "glue ball" chiral 
superfield 

S = -tilWWa] = tr[A ■ A] + ■ ■ • . (2.4) 

The Kahler potential which determines the kinetic energy terms of the fields T, B, and B 
was hypothesized in the work of Masiero and Veneziano [?, Q to take the following form 

K[T,B,B] = ArtilT^T] + Ab{B^ B + B^B). (2.5) 

Our main results will rely on weaker assumptions about the Kahler potential, in particular, 
that it is nonsingular on the space of supersymmetric vacuum states. However, we will 
support our general remarks by explicit calculations using this simple model. We expect 
(^.51) to be the correct form of the Kahler potential near the origin of moduli space, in the 



cases for which the mesons and baryons give an effective infrared description of the theory. 
This is not the case for Nc > Nf^. 

2.1 Soft Super symmetry Breaking 

We break supersymmetry by adding mass terms for the squarks and gaugino fields, 

AC = -mKlQf + \Qf) + {m,S + h.c), (2.6) 



where, in ( p.6|) , Q, Q, and S are the scalar component fields of the superfields. The squark 
mass term is the unique soft supersymmetry breaking term which does not break any of the 
global symmetries (|2.1| ) of the original model. The gaugino mass term breaks only the U {1)r 
symmetry, and thus breaks the global symmetry of the supersymmetric model down to that 
of ordinary QCD with Nf massless fiavors. Actually, there are two consistent possibilities 
for SSB theories which we denote by R and ^ 

R : rriQ 7^ rUg = =^rn'^ -»oo standard QCD with Nf quarks + adjoint fermions 



I? : niQ 7^ rrig 7^ ^^m^ standard QCD with Nf quarks 

Since we will be working in the language of the low-energy effective Lagrangian, we must 
ask how the supersymmetry breaking term ( |2.(j| ) shows up in this Lagrangian. To work this 
out, rewrite ( |2.6| ) in the superfield form 

A£ = y d^eMQ{Q^e^Q + Q^e~^^Q) + J d^OMgS + h.c. (2.7) 

where Mq is a vector superfield whose D component equals {—mQ) and Mg is a chiral 
superfield whose F component equals mg. It is straightforward to see that these superfields 
are gauge-invariant and neutral under all of the global symmetries. 

The effective Lagrangian description of AC for Nf < A^^c + 1 is then given by writing 
the most general Lagrangian built from T, B, B and a fixed number of factors of Mq and 



Mg. The supersymmetry breaking terms have an ambiguity related to that of the Kahler 
potential, because many possible invariant structures can be built from T B, and B. In 
our explicit calculations, we will assume that the coefficient of Mq is quadratic in these 
fields. This assumption holds near the origin of moduli space. Then the first order soft 
supersymmetry breaking terms in the effective Lagrangian are 

AC = Jd^e^BrMQiilT^T] + BbMq{B^B + B^B} + M'gM{T,B,B) + h.c.) 

+ J d'^eMg{S) + h.c, (2.8) 

clarifications and remarks: 

(i) A4{T, B, B) is a function of the effective Lagrangian superfields which is neutral 
under the global symmetries. 

(ii) The quantity (S) in (|2.8|) should be a combination of the effective Lagrangian chiral 
superfields which has the quantum numbers of S. In general, this condition restricts that 
function to be proportional to the expectation value of S as determined from the effective 
Lagrangian of refs. 0, 0| which includes S" as a basic field. 

(iii) The squark mass terms in (|2.8| ) are not the most general terms that can be written 
down. Higher order terms in the fields, suppressed by powers of A, may appear. However, 
we expect ( p^ ) to be approximately true near the origin of moduli space T = B = B = 0. 

(iv) We have assumed throughout this work that the coefficients Bb and Bj- are pos- 
itive. This is not an innocent assumption since had they been negative chiral symmetry 
breaking would have taken place instead of several of the exotic features. At present we do 
not have a proof of the positivity of those coefficients, but the picture we derive assuming 
that they are positive seems to be consistent. 

(v) The ratio of coefficients Bb/ -Bt will be important to our later analysis, but this ratio 
cannot be determined from the effective Lagrangian viewpoint. At best, we can argue naively 
that the coefficient of the mass term of a composite field should be roughly proportional to 
the sum of the coefficients of the mass terms of the constituents. This would give the relation 
Bb ~ which the reader might take as qualitative guidance. 

(vi) It is important for our analysis that the behavior of the theory is non-singular 
when adding the squark and gluino masses, i.e. that no new non-perturbative effects occur. 
In general it is not possible to prove this in non-supersymmetric theories, but a proof of this 
is possible in softly broken supersymmetric theories, when the soft breaking can be viewed 
as arising via spontaneous breaking of supersymmetry. [^] 

3 Baryon number violation (Nf = Nc) 

The cases of Nc > Nf are discussed in ref. here we start with the case Nf = N^. The 
low-energy effective Lagrangian of the supersymmetric limit contains both meson and baryon 
superfields. The quantum numbers of the T*j, B and B are summarized in the following 
table. 
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This model has a manifold of quantum mechanical supersymmetric ground states, in 
which the meson and baryon fields satisfy the relation (in units where A = 1) 



detT - BB = 1 . 



(3.9) 



Several forms for the superpotential are consistent with this relation. The S'-dependent 
superpotential, for example, has the form W = 5'log(detT — BB) . Note that this super- 
potential leads to conditions for a supersymmetric vacuum state which imply not only (|3.9| ) 
but also the constraint (S) = 0, so that the U{1)r symmetry is not spontaneously broken. 

The presence of a manifold of degenerate vacuum states not related by a global sym- 
metry is necessarily accidental unless it is a consequence of supersymmetry. Thus, any such 
degeneracy should be broken as soon as supersymmetry breaking terms are added to the 
Lagrangian. To first order, this is the main effect of the soft supersymmetry breaking per- 
turbation. To analyze this effect, we should restrict our attention to the values of T, B, and 
B obeying the constraint ( |3.9| ) , for which the vacuum energy vanishes in the supersymmetric 
limit, and study the behavior of the supersymmetry breaking potential over this space. 

For the R models, the soft supersymmetry breaking terms ( |2.8|) lead to the potential 



AV" = BTmltT[T^T] + BBm^iB^B 



B^B) 



(3.10) 



Using SU{Nf) X SU{Nf), we can diagonalize T to complex eigenvalues tj. Parameterize the 
baryon fields as B = xb , B = —^b with x and b complex. This potential has three types of 

stationary points 

^ 4 \ 
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with X = 1 in all cases. The method of effective Lagrangians cannot tell us which of the two 
vacuum states at 6 = and T = is the preferred one. This depends on the ratio Bb/Bt, 
which is a phenomenological input to the effective Lagrangian analysis. The vacuum at 
T = is locally stable if Bb > Bt and is globally stable if Bb > {Nf /2)Bt. Our naive 
estimate puts the theory just at the boundary at which the two vacuum states have equal 
energy. 

The structure of the spectrum in the 6 = and T = cases is summarized as follows 
(for details see 0). 
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Notice that since ( ^.91 ) is a superfield constraint, it also removes one fermion from the 
theory, specifically, the fermionic partner of det[T] or of B + B. 

For the ^ models with small rrig the qualitative picture remains the same. We showed 
earlier that the superpotential implies that, in the manifold of supersymmetric vacuum states 
about which we are perturbing, (5*) = 0. Thus, the superpotential term proportional to Mg 
does not contribute to the vacuum energy. More generally, since Mg, T, B, and B are all 
invariant under while a superpotential has R charge 2, this term does not contribute 

to the superpotential to any order in mg. There are possible Kahler potential terms involving 
Mg. However, near the vacuum with 6 = 0, these will be polynomials in B and B of order 
at least 2, and near the vacuum with T = they will be polynomials in T of order at least 
2. Thus, these terms will not affect the presence of stationary points of the vacuum energy 
at these positions in the field space. Examples of D-term terms that can be induced are 



AC = / d^9Mg{CT det T + CbBB) + h.c. 



(3.11) 



where Ct and Cb are some constants. If one begins from the effective Lagrangian including 
S, with the canonical superpotential and Kahler terms. 



C = J d^eS*S + J d^e{S logidetT - BB) + MgS) + h.c. 
and integrates out S, one finds a breaking term 

AC= I d^eWglogidetT - BB) + h.c. 



(3.12) 



(3.13) 



which gives qualitatively similar results. In the following discussion, we will work with p.ll| ). 
Once the f/(l)R is broken it does not protect fermions from acquring mass. Around 6 = 
the fiavor adjoint and baryonic fermions, acquire the following masses m^j = \^mg and 

m^B = cind no zero-mass fermions remain. Around the T = vacuum ipTj and ifjr 

remain massless but the baryonic fermion aquires a mass m^B 



Ab "^9- 



3.1 Toward the Decoupling Limit 



The transition from the region of weak supersymmetry breaking to the decouphng hmit 
mg, nig — > oo is very different in the four models of R and ^ with vacua at 6 = and T = 0. 

(i) ft b=0- In this case (hke in the cases of Nc > Nf) the global group SU{Nf) x 
SU{Nf) X U{1)b is broken spontaneously to SU{Nf) x U{1)b, leaving no massless parti- 
cles except for the required Goldstone bosons. It is reasonable to expect that there is a 
smooth transition from the situation of weak supersymmetry breaking to the decoupling 
limit rriQ, rrtg oo. In QCD, chiral symmetry breaking is characterized by a nonzero vac- 
uum expectation value of the quark-antiquark bilinear, i^qi^j in our present notation. In 
the language of the supersymmetric effective Lagrangian, this operator is a part of the F 
term of the superfield T^j. The expectation value of this term may easily be found to 

Nc 

be proportional to rriQ " ^ . Thus, the F term of T does obtain an expectation value in 
the vacuum state that we have found. This expectation value naturally becomes a nonzero 
expectation value for the quark bilinear in the decoupling limit. As mq increases, the quark 
bilinear becomes larger while the squark bilinear becomes smaller, in exact accord with our 
expectations. 

(ii) R b=0- As itlq is taken to infinity, the squarks decouple, and the model becomes 
a purely fermionic QCD theory with Nj quark flavors plus one fermion flavor in the adjoint 
representation of the gauge group. For small values of the supersymmetry breaking mass 
ttIq, this vacuum contains massless fermions corresponding to the fermionic components of 
the superflelds T, B, and B. We might think of these as being built out of scalars, with one 
squark replaced by a quark to give the composite spin i. But it is also possible to build 
composite bound state with the same quantum numbers purely out of fermions, by replacing 

Q' - ^"-^a , Qj - A>Q„,. , (3.14) 

where a is a two-component spinor index and the gauge indices are implicit. Notice that 
this combination has the same quantum numbers as the squark, including zero R charge. 
Then, for example, the fermion created by T^j could be constructed as 

i^Ta'j - i^hc^^^p, ■ (3-15) 

With this replacement, the composite fermions are built only out of constituents which 
remain massless as the squarks are decoupled. Thus, it is apriori reasonable that the 6 = 
vacuum of the R model could go smoothly into a vacuum of the purely fermionic QCD theory 
described above. This vacuum would have broken SU{Nf) x SU{Nf) but unbroken chiral 
zero values for the vacuum expectation values of quark-antiquark bilinears, massless 
composite fermions in the adjoint representation of flavor SU{Nf), and massless baryons. 
We will refer to this scenario as 'option 1'. 

The other possibihty for this model is that, after the squarks decouple, the gluino flelds 
pair-condense, in a second-order phase transition at some value of m^, and the nonzero value 



of the condensate (A • A) spontaneously breaks f/(l)j?. In this case, the physics would revert 
to the usual symmetry-breaking pattern of QCD, and the composite fermions would become 
massive. The gluino condensate would make itself felt only by providing an extra SU{Nf)- 
singlet Goldstone boson. We will refer to this scenario as 'option 2'. 

(iii) /? T=0- In this model, in the supersymmetric limit, the massless composite 
fermions belong to the {Nf,Nf) representation of an unbroken flavor group SU{Nf) x 
SU{Nf). There are no constraints from the symmetry, which is explicitly broken, 
or from baryon number, which is spontaneously broken. With this freedom, can we build 
these fermionic composites out of fields that survive in the decoupling limit mQ,mg oo? 
For Nc even, it is impossible, because the only constituents available are the quarks ipq, i'gj, 
and gauge-invariant states must contain an even number of these. For Nc odd, however, it 
is possible to build composites with the correct quantum numbers, as follows: 

^PTa', - e^'-' 6,,...,V^g ■ ■ ■ V^^^, , (3.16) 

where -ipQ is the right-handed fermion field in (Q)*. The {Nc — 1) right-handed fermion fields 
must be contracted into a Lorentz scalar combination. For the case Nf = Nc = 3, eight of 
the nine fermions in ( p.l6| ) have the quantum numbers of the baryon octet in QCD. 

However, in this case, there are two compelling arguments that the spectrum which we 
find cannot survive to the decoupling limit. In the limit rriQ oo, even without introducing 
nig, we have a vectorlike gauge theory of fermions. For such theories, the QCD inequalities of 
Weingarten [|lOl and Vafa and Witten [jlll apply. In we use Weingarten's method to prove 
that, in the decoupling limit, flavor nonsinglet composite fermions must be heavier than the 
pions, which are massive in the T = vacuum. Alternatively, we can apply the theorem 
of Vafa and Witten in the decoupling limit to show that vectorlike global symmetries, in 
particular, baryon number, cannot be spontaneously broken. 

By either argument, the T = vacuum state must disappear via a second order phase 
transition at a flnite value of rriQ. Most likely, this vacuum becomes locally unstable with 
respect to a decrease in the expectation value of 6, driving the theory back to the more 
familiar vacuum at 6 = 0. 

(iv) R T=0- In the supersymmetric limit we have massless fermions in the following 
representations {Nf, Nf, —1) -|- (1, 1, —1) of the unbroken symmetry group, corresponding to 
the fermions in T and a linear combination of the fermions in B and B. Both multiplets 
are necessary to satisfy the anomaly conditions involving U{1)r. The arguments that we 
have just presented for the T = vacuum in the ^ models apply equally well to the R case. 
Again, we must have a second-order transition, probably with an instability to the 6 = 
vacuum. There are then two possible endpoints, depending on which option is chosen for 
the 6 = vacuum. If the option 1 for the 6 = vacuum is correct, it is not necessary that 

be spontaneously broken in this transition. 



4 No x^Si^, massless fermions (Nf = Nc-\- 1) 

In the case of Nf = N^ + l, the low-energy effective Lagrangian in the supersymmetric hmit is 
expressed in terms of the baryon, anti-baryon and meson superfields which transform under 
the global symmetry as follows 



B : (iV;,l)i,i-^ B : il,Nf)_,^,_^ T : {Nf,Nf\^ (4.17) 

where the second subscript is the R charge of the scalar component of the superfield. In the 
supersymmetric theory the low energy effective theory is described (at least near the origin 
of moduli space) by the Kahler potential given by ( p.5| ), and by the following superpotential 

W = BiT]B^ -detT. (4.18) 
The supersymmetric vacuum is, thus, described by a moduli space characterized by 



B{r; =0 T^B'= ±-e,^_,^/^-^-fTil...T^^ - B.^ B'^f = 0. (4.19) 



As was argued in |T^ , these equations correctly describe the moduli space of vacuum states 
in the full quantum theory. At the origin of the moduli space, < T >=< B >=< B >= 0, 
where the full global symmetry ( p.l|) remains unbroken, there is a further consistency check 
for the low energy behavior. The fermionic components of the low-energy superfields ( [4.17] ) 
match the global anomalies of the underlying theory. 

When we break supersymmetry by squark and gluino masses, we add to the effective 
Lagrangian the mass terms for T, B and B indicated in ( p.8|) . Since we are adding terms to 
the potential which are positive and vanish at the origin of moduli space, it is obvious that 
the origin becomes the only vacuum state of the theory. All of the scalar particles in the 
effective theory obtain mass terms proportional to Btttlq or Bbttlq. 

Though all of the scalars obtain mass, all of the fermions remain massless. The su- 
perpotential (|4.18|) is a least cubic in the fields, thus, any mass term derived from this 
superpotential vanishes at the origin. Similarly, in the ^ case, the Mg term in ( p.8| ) requires 
a function of T, B, and B which is neutral with respect to the global group; the only such 
functions, quadratic in fields, are tr[T"l'T], B'^B, and B^B, and these do not give rise to 
fermion masses when multiplied with Mg. In fact, it is required that no fermions should ob- 
tain mass, since the full multiplet of fermions in T, B, and B is needed to satisfy the 't Hooft 
anomaly conditions for the global symmetry group SU{Nf) x SU{Nf) x U{1)b x U{1)r. 



5 Seiberg's duality in SSB models (Nf >Nc + 2) 

No solution of the 't Hooft anomaly matching conditions for SQCD involving gauge invariant 
bound states is known for Nf > {N^ + 1). However, Seiberg has suggested a compelling 



solution to these constraints in terms of a novel type of duality in which new gauge degrees 
of freedom are dual to the original quarks and gluons |l|]. The pair of dual theories are 
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On top of the fact that these theories have the same 't Hooft anomalies, they were shown to 
be equivalent in the infrared also by analyzing their "chiral rings" , fiat directions and mass 
perturbations. In the supersymmetric limit of the magnetic theory the scalar potential has 
a moduli space of vacua, which includes the point < T >=< q >=< q >= at which the 
chiral symmetry is unbroken 

Now add squark masses to the theory. Their effect should be seen in the effective 
Lagrangian. In case that the latter is given by the dual theory we should add to the effective 
Lagrangian of the dual theory the term 

AV = BrmliiiT^T) + 5,mJ(|gp + Igp), (5.20) 

at least near the origin of moduli space. After we add this perturbation, the only minimum of 
the potential is at < T >=< q >=<q>= 0. Thus, adding a squark mass leaves the theory 
in the phase in which the chiral symmetry is unbroken. All scalars get masses (originating 
only from AV^, since the original scalar potential is quartic in the fields), while all fermions 
remain massless. As in the original supersymmetric theory, this complement of massless 
fermions has just the right quantum numbers to satisfy the 't Hooft anomaly conditions for 
completely unbroken chiral symmetry. 

The glueball operator tr(lV^) is identified (up to a sign) between the original and the 
dual theory [|l^]. Thus, to leading order in rrig, a gluino mass in the original theory is just 
equal to a gluino mass in the dual theory. Adding this term breaks the ?7(1)_r symmetry, 
but the SU{Nf) x SU{Nf) global symmetry still remains and protects the dual quarks from 
getting a mass. Thus, we find the same spectrum in the R and ^ cases, except that in the 
latter case the dual gluino, which can be an asymptotic particle, becomes massive. 

Let us discuss now the infrared description of the theory. For small rriQ (and small rrig, 
in the ^ case) the situation is schematically described in the following diagram 
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The picture describes a scenario where in the range between the : signs, namely < 
Nf < for the R case and between the * signs, namely, ^N^ < Nf < ^N^ for the ^ case, 
the two theories flow to a common fixed point in the IR. In the range Nf < ^-Nc {Nf < ^N^) 
the R (^) system flows to the free magnetic theory in the IR, whereas for Nf > {Nf > 
^Nc) the picture in the infrared for R (^) is that of the free electric theory. 

The discussion of the decoupling limit for these theories is very similar to that for the 
Nf = Nc theory. Thus, for those values of Nf which have massless fermions for small values 
of rriQ we must have a second-order phase transition as rriQ is increased. It is not clear 
how the theory behaves on the other side of this phase transition. In the |? case obviously 
only a SU{Nf) x U{1)b symmetry remains, with no massless fermions. This situation is 
schematically described as follows 
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6 Problems of Approximate Supersymmetry on the 
Lattice 



(i) Can the phenomena we have discussed in this work be seen in lattice gauge theory 
simulations? In general, gauge theories on the lattice cannot be made supersymmetric at 
the fundamental level. We expect that lattice simulations of these theories will also contain 
small dimension 4 perturbations which violate supersymmetry. Our analysis has been based 
on the assumption that, if the phenomena discussed by Seiberg survive perturbations which 
are relevant in the infrared, they should also survive small marginal perturbations. 

(ii) Can one reach the continuum limit in such lattice simulations? Typically in lattice 
gauge theory simulations with scalar flelds, there is no continuum limit; instead, one flnds a 
flrst order phase transition as a function of the scalar fleld mass parameter WM. Renormal- 



ization effects in a gauge theory can induce an unstable potential for a scalar fleld coupled 
to the gauge bosons, leading to a 'fluctuation- induced flrst-order phase transition'. How can 
this be avoided? 



(iii) To analyze this question, consider the renormahzation group equations for an 
approximately supersymmetric gauge theory. Viewed as a conventional renormalizable gauge 
theory, SQCD has three coupling constants, the gauge coupling g, the quark-squark-gluino 
coupling and the four-scalar coupling go- The scalar potential has the specific form 
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(6.21) 



where r"^ is an SU{Nc) matrix. If we relax the constraint of supersymmetry, there are four 
possible invariants under the symmetries of the problem, including the continuous global 
symmetries and parity Q ^ q\ The most general linear combination of these invariants 
can be generated by the renormahzation group flow. 

We can view the effects of renormahzation relatively simply, however, by restricting 
our attention to the vicinity of the surface given by the three couplings g, gx, and go- In 
this surface, the beta functions of the three couplings are given by 



1 



(47r)2 
1 

1 

(4^ 
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gx9\m + 3C2(iV,)) - gl{3C2{N,) + Nf 
4gtN, + 2gUK - Nf - 2C2(iV,)) + Ug^g^C^m - SglgfC^m 



where C2{Nc) = (NJ^ — l)/2Nc. These three functions all reduce to the standard SQCD 
beta function on the supersymmetric subspace; for g^ = g\ = g^, Pg = Pg^ = (3g2^/2g. 
Note that, for ~ Nf and the three couplings in reasonable ratio, all three couplings are 
asymptotically free. 

The potential instability to a first order phase transition arises because a new structure 
in the potential is induced by the renormahzation group flow. To lowest order, the form of 
the potential induced is 



9l 



On the surface g% = ^1 the beta function for g'\, is 

1 



A 



9l 



(47r) 



g\> 



(6.22) 



(6.23) 



This equation implies that, if one leaves out the gluinos, g\ becomes negative in the infrared, 
leading to a fluctuation-induced first-order phase transition. According to ( |6.23| ), this effect 
is removed if the lattice simulation includes gluinos, and if the gluino coupling g\ is large 
enough. 



With this provision to avoid possible first-order phase transitions, we expect that lattice 
simulations with an approximately supersymmetric action can reach the continuum limit and 
test our predictions for softly broken supersymmetric QCD. 
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